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Abstract 

We study the asymptotic solutions of the Schrodinger equation for the color-singlet 
reggeon compound states in multi-color QCD. We show that in the leading order of asymp- 
totic expansion, quasiclassical reggeon trajectories have a form of the soliton waves propa- 
gating on the 2-dimensional plane of transverse coordinates. Applying the methods of the 
finite-gap theory we construct their explicit form in terms of Riemann theta-functions and 
examine their properties. 
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1. Introduction 



It widely believed that the Regge asymptotics of hadronic scattering amplitudes in high-energy 
QCD should be described by an effective Regge theory. In this effective theory reggeized gluons, 
or reggeons, play the role of a new collective degrees of freedom. Reggeons form a color singlet 
compound states, QCD Pomerons and Odderons, which propagate in the t— channel between 
scattering hadrons and contribute to a power rise with energy of the physical cross-sections IJ. 
In the Bartels-Kwiecinski-Praszalowicz approach (BKP) 0, the color singlet reggeon compound 
states are built from a conserved number N = 2,3, ... of reggeons. For a given number of reggeons, 
N, the wave function of these states xn depends only on the transverse reggeon coordinates 
b±,j = (xj,yj) (j = 1, ...,N) and it satisfies the BKP equation 

T-Ln\xn) = £n\xn) (1.1) 

where Ti^ is an effective QCD Hamiltonian acting on the transverse coordinates of N interacting 
reggeons. QCD Pomeron and Odderon appear as the states in the spectrum of Ti^ with the 
maximal energy en- 

The BKP equation has a number of remarkable properties in multi-color limit of QCD, N c — > 
oo and a s N c = fixed. Firstly, introducing holomorphic and antiholomorphic coordinates on 2-dim 
plane of transverse reggeon coordinates, Zj = Xj + iyj and Zj = Xj —iyj, respectively, one can find 
that the wave function xn splits into the product of holomorphic and antiholomorphic parts 

Xn{xi,Vi, ...,x n , y N ; x , y ) = (f N {zi - z , z N - z ) ■ Tp N {z x - z , z N - z ) , 

where (xo,yo) is the coordinate of the center-of-mass of the N reggeon compound state. The 
effective hamiltonian T-Ln is invariant under SL2 transformations 

^azj + b ^ azj + b 

Zj ~* czj + d 1 Zj ^ czj + d' ^ ' ' 

with ad — be = ad — be = 1 and j = 0, 1, N, while the wave function is transformed as 

Xn -> (cz + d) 2h {cz + df h XN . (1.3) 

The conformal weights take the values h = ly 2 — iv and h = 1^ — iv with m = 7L, v = 1R 
corresponding to the principal value representation of the SL(2,€) group. Secondly, the QCD 
effective Hamiltonian Ti^ is closely related to the Hamiltonian of XXX Heisenberg magnet of 
noncompact SL(2,€) spin s = and, as a consequence, the system of iV interacting reggeons 
possesses a large enough family of conserved charges in order for the Schrodinger equation 
( |1 . 1|) to be completely integrable 0, ||. This implies that instead of solving the BKP equation 
one can define the holomorphic wave function tp^ as a common eigenfunction of the family of 
N — 1 mutually commuting operators q n ] = 

Qk = z hh z hh--- z ikh r Ph r Ph---'Pi k 1 k = 2,...,N (1-4) 

i<ii<i2<...<jfc<A r 

with z a b = z a — z\ } and pj = —i-^r- There are also additional constraints on (p^ 

(L 3 + h)ip N (z 1 ,...,z N ) = 0, L+<p N (zi, ...,z N ) = 0, (1.5) 
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which follow from the transformation properties of the wave function (1.3). Here, the notation 
was introduced for the SX 2 generators 



N N N 

k=l k=l k=l 



with pk = —idk- The conserved charges qu defined in ( |1.4D have the meaning of higher Casimir 
operators of the SL(2,<C) group those eigenvalues determine the quantum numbers of the N 
reggeon compound states. In particular, 

Q2 = - z %PjPk = -L3L3 - \(L_L + + L+LJ) 
j>k A 

is a quadratic Casimir operator and its eigenvalue corresponding to the wave function (p^ is given 
by qi = —h(h — 1) provided that conditions ( |1.5| ) are fulfilled. 

To find the holomorphic reggeon wave function (p^ one has to diagonalize simultaneously the 



remaining N — 2 commuting hamiltonians qs, q^ on the space of functions satisfying ( |1.4|) . 
This leads to a system of N coupled partial differential equations on the wave function (p^. The 
same system can be also interpreted |J as a set of Schrodinger equations, in which the conformal 
weight h plays the role of the effective Planck constant, h = 1/h. This suggests to consider the 
quasiclassical limit of large h and develop the WKB expansion for the holomorphic reggeon wave 
function 

<p N (zi, z N ) = exp(iS + iS 1 + ...) , (1.6) 

where Sq = 0(h), Si = O(h ), ... at large h and Sk = Sk(zi, ...,2jv) = 0{h}~ k ). In a complete 
analogy with quantum mechanics, the leading term, Sq, defines the action function for the system 
of N reggeons. For a given set of quantum numbers q^, it satisfies the system of the Hamilton- 
Jacobi equations 



J2 z k^= ih , Zj 2: fc^ = ' Qa[z,^-) = q a , a = 3,...,N, (1.7) 




where z = (z\, z N ) and q a {z,p) stands for the symbol of the operator ( |1.4|) . 

For N = 2 reggeon state, the BFKL Pomeron, the Hamilton- Jacobi equations ( |1 . 7| ) can be 
easily solved 

S (z 1 ,z 2 ) = -ih\n (---), (p 2 (z uZ2 )= (—) H x h + Oih- 1 ) 

\Zi Z 2 J \Z\ZiJ L 

The last expression defines the holomorphic wave function in the leading order of the WKB 
expansion. Comparing it with the well-known expression for the wave function of the BFKL 
Pomeron we find that it is exact to all orders in 1/h. 

For N > 3 reggeon states, the Hamilton- Jacobi equations ( |L7| ) can be solved in the separated 
variables || . Their solutions describe the compound N reggeon states as a system of TV classical 
particles moving along quasiperiodic trajectories. The corresponding action-angle variables were 
constructed in 0, where the close relation between classical reggeon trajectories and finite-gap 
soliton solutions was established. The main goal of the present paper is to apply algebro- 
geometrical methods [|1(J in order to obtain the explicit form of the reggeon trajectories on the 



2— dimensional plane of transverse coordinates in terms of Riemann theta-functions JTT 
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2. Hamiltonian flows 

In the leading order of the WKB expansion, we replace in ( |1.4|) the momenta operators p k by the 
corresponding classical functions to get the family of N — 1 classical Hamiltonians q%, ?j\r- For 
the system of N reggeons with the coordinates z k , momenta pk and the only nontrivial Poisson 
bracket {xk,p n } = 8kn, these hamiltonians generate the hierarchy of the evolution equations 

^ = f z Q } = ^ ^ = h a } = (2 1) 

dt a {Zn,q ^ dpj dt a iPn ' q ^ dzj 1 j 

with t a (a = 2, N) being the corresponding evolution times and d ta q n = {q n , q a } = 0. Their 
solutions define the reggeon trajectories z k = z k (t 2 , ...,tv) subject to the periodicity condition 

z k+N (t) = z k (t) , Pk+N{t) = Pk{t) ■ (2.2) 



The system of evolution equations ( |2.1|) has a sufficient number of the integrals of motion q k to 
be completely integrable. It admits the Lax pair representation, which can be found using the 
equivalence of the system of N reggeons and the XXX Heisenberg magnet of SL 2 spin 0. Namely, 
for each reggeon we define the Lax operator as 



with E being an arbitrary complex spectral parameter. Then, the evolution equations (2J) are 
equivalent to the matrix relation 

d ta L k = {L k , q a } = A k %L k - L k A k a) , (2.4) 

where A k °^ (E) is a 2 x 2 matrix depending on the choice of the hamiltonian. As an example, for 
the hamiltonian q^ one can obtain 

A[ N) = E^L- ( 1 U (^, -1) (2.5) 



Z k -l,k \ z k 

and the corresponding evolution equations look like 

a Q N a ( l 

Ot N x k = — , d tN p k = q N 



Pk \Z k -i yk Z kjk+ i ) 

The explicit form of A k for the remaining hamiltonians can be deduced from the Yang-Baxter 
equation for the Lax operator L k and it will not be important in the sequel. 



The integration of the evolution equations ( |2.1| ) and (|2.4j ) is based on the Baker-Akhiezer 



function ^ k (E;t 2 , •••,£tv)- It is defined as a solution of the following system of matrix relations 

L k (E)V k = V k+l , d ta V k = A k a \E)y k , tf fc =^M, (2.6) 
where <fi k and Xk are scalar components. We construct the monodromy matrix 

r(B ) = L w ( B )... £ ,( B)£l ( B )=gg52) (2.7) 



and verify using ( |2.4|) that its eigenvalues are integrals of motion due to dt a T(E) = Af T—TA " . 
The monodromy matrix generates the shift T(E)^>i = ^n+i and we impose the periodicity 
condition on the Baker- Akhiezer functions by requiring ^/^(E) to be the Bloch-Floquet function 

y k+N (E)=e P{E) *k(E). (2.8) 

Here, P{E) is an eigenvalue of the monodromy matrix ( |2.7| ) and it satisfies the characteristic 
equation 

det(T(£) - e p ^) = e 2P ^ —A(E) e p ^ +1 = 0, (2.9) 

where A(E) = tr T(E) can be calculated using ( |2.7p and fl2.3| ) to be a polynomial of degree A" in 
the spectral parameter with the coefficients given by the integrals of motion 

A(E) = 2 + q 2 E 2 + q 3 E 3 + ... + q N E N . 

Introducing the complex function y(E) = (e p ^ — e~ p( - E ^)/E one rewrites ( |2.9| ) in the form of 
algebraic complex curve [] 

T N : y 2 = E- 2 (A 2 (E)-4) = (q 2 + q 3 E + ... + q N E N - 2 )(4 + q 2 E 2 + q 3 E 3 + ... + q N E N ) . (2.10) 

For any complex E in general position the equation ( |2.9| ) has two solutions for P(E), or equiva- 
lently y(E) in (|2.10|) . Under appropriate boundary conditions on ^k(E) (to be discussed later) 
each of them defines a branch of the Baker- Akhiezer function, ty^(E). Then, being a double- 
valued function on the complex E— plane, ^j^(E) becomes a single- valued function on the Rie- 
mann surface corresponding to the complex curve T^. This surface is constructed by gluing 
together two copies of the complex E— plane along the cuts [e2,es], [e2jv-2>ei] running be- 
tween the branching points ej of the curve (see for detail Fig. 1), defined as simple roots of the 
equation 

A 2 (e,)=4, j = l,...,2iV-2. (2.11) 

Their positions on the complex plane depend on the quantum numbers q 2 , q 3 , qn of the 
reggeon compound state. In general, thus defined Riemann surface has a genus g = N — 2 and it 
depends on the number of reggeons, e.g. it is a sphere for N = 2 state, or BFKL Pomeron, and 
a torus for A^ = 3 state known as QCD Odderon. 

To distinguish points Q = Q(E) on the Riemann surface lying above E on the upper and the 
lower sheets we assign them a sign Q = (±,E), where the upper (+) sheet corresponds to the 
asymptotics y ~ q^E 1 ^' 1 as E — » oo. In particular, one finds from Q2.9|) the behaviour of the 
Bloch multiplier on different sheets of Tn for E — > oo as 

e ±p ( £ ) = q N E N x (1 + 0(1/ E)) , (2.12) 

where ± sign corresponds to the choice of the sheet. 

2.1. Boundary conditions 

Being written in components, the first condition ( |2.6| ) on the Baker-Akhiezer function looks like 

Xn+l - Xn = Z n {(f n+1 - ip n ) = Ez n p n (x n ~ Z n <f n ) , (2.13) 

1 To make a correspondence with notations of ||, one has to redefine the local parameter on the curve as 
E=l/x. 
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Figure 1: The definition of the canonical basis of cycles on the Riemann surface IV. The dotted 
line represents the part of the f3— cycles belonging to the lower sheet. Cross denotes a projection 
of the points Qq G on the complex plane. The path 7q-q+ does not cross the canonical 
cycles and goes from the point Qq on the lower sheet to the point Qq on the upper sheet. 



where (p n and \ n are single- valued function of Q(E) on the Riemann surface T N . We notice that 
for E = the solutions to this equation do not depend on the reggeon number n. Moreover, 
using ( |2.3| ) and (|2.5| ) one can show that L k (E = 0) = 1, A£*\E = 0) = and, therefore, there 



exist two special points on the Riemann surface Q^ = (±, E = 0), at which the Baker- Akhiezer 
function takes constant values, ^k(Qo', {t}) = const. The choice of the constants implies the 
normalization of the Baker- Akhiezer function. 

Let us show that the values ^(Qo) are fixed by the constraints ( |1.5|) . At the vicinity of 
E = on one the sheets of the Baker- Akhiezer function can be expanded as 

* n (Q; {*}) = ¥<°> + E + 0(E 2 ) , 

where coefficients depend on the choice of sheet and the leading term \l/( ) is a constant. Let us 
substitute this identity into relation ^i + n = T(E)^i = exp(P(E))^/ 1 and expand its both sides 
in powers of E. Taking into account small-E expansion of the monodromy matrix 



T(E) = l + iE\ L T Z L -)+0(E 2 ), (2.14) 



one obtains that is an eigenvector of the following matrix 




where conditions ( |1.5| ) were imposed. We find two solutions for corresponding to P'(0) = 
±ih, which fix (up to overall constant) the value of the Baker- Akhiezer functions at the points 
Qq as 

*n(Q + ) = (iV *n(Qo)=(^PV (2-15) 



0, 

N 



Here, h is the conformal weight of holomorphic wave function and L_ = —iJ2k=iPk- in th e 
quantum case, the operator L_ commutes with the hamiltonians q^, but it does not take any 
definite value on the subspace (|L~5|), since [L_,L 3 ] = 2iL_ and = iL 3 . In the quasiclas- 

sical limit, L_ has a zero Poisson bracket with and it does not depend on the evolution times 
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tk- Its value is fixed by the initial conditions as 



N 



£ Pfc (0) = -zP t 



tot ) 



k=l 



with P tot being the total holomorphic momentum of N reggeons.0 

We notice, however, that the relations (|2.6| ) and ( 2.13|) are invariant under SL 2 transforma- 
tions 



Pk -> Pk(cz k + df , 



d c 
b a 



with the coordinates z& transformed as in ( |1.2|) . Since this transformation changes the value of 
L_, one can choose its parameters in such a way that one puts L_ = and keeps the relations 
L 3 = —h and L + = unchanged.^ The relation between original and transformed reggeon 
coordinates looks as follows 



—l- 



(0) 



2h 



(0) 1 
Zn 



Pn = Pn 



p 

-, •-'tot (0) 



(2.16) 



,«)) 



It allows to integrate the evolution equations (|2.1|) for L_ = to get expressions for z n and p n 
and then restore the physical solutions using (|2.16|) . 

Therefore, in what follows we put L_ = and obtain the normalization conditions ( [2.15| ) for 
the Baker-Akhiezer function in the canonical form (up to overall constant factor) 



(2.17) 



It easy to see that for L_ = L + = the solutions of the evolution equations (|2.1| ) and (|1.5| ) are 
invariant under rescaling of the coordinates 



(0) 
Zn 



(0) 
pZn 



P 



(0) 



-1 (0) 
P Pn 



(2.18) 



For the Baker-Akhiezer function, this corresponds to a freedom in choosing an overall normal- 
ization factors in ( |2.17| ). One can further simplify analysis of the evolution equations ( [2.1| ) for 
zh 0) and p { n by noticing that qi = —L 3 L 3 = —h 2 and the dependence on the "lowest" time t 2 is 
given by d t2 zh 0) = — 2L 3 |^- = 2ihz n °\ and similarly for p { n\ leading to 



Zn 0) (t 2 ,t 3 ,...,t N ) =Zn°\0,t 3 ,...,t N )e 



(0), 



^2iht 2 



p l n\t 2 ,t 3 ,...,t N ) = P { n , (0,t 3 ,...,t N )e 



-2iht 2 



(2.19) 

Thus, one can forget for a moment about t 2 — dependence of the reggeon coordinates and restore 
it in the final expressions using ( |2.19| ). 

2 In hadronic scattering amplitudes, Ptot is defined as a holomorphic component of the total momentum trans- 
ferred in the t— channel. 

3 These three conditions form a set of the second-class constraints for the TV reggeon system. Their quantization 
leads to the Baxter equation for the wave function of the N reggeon compound state |0| . 
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2.2. Singularities of the Baker- Akhiezer function 

Let us show that the components of the Baker- Akhiezer function ^> n {Q) are meromorphic func- 
tions on the Riemann surface Tjy having g + 1 = N — 1 poles, the same number of zeros and 
essential singularities at two infinities = (±, oo) situated on the upper and lower sheets of 
Tjv. Let us take the B— component of the monodromy matrix ( |2.7p and observe that B(E)/E is 
a polynomial of degree N — 1 in the spectral parameter E and B(E) = iL^ E + 0(E 2 ) as E — > 
according to (|2.14|) . Then, one defines the points E\, -Etv_i as roots of this polynomial 



B(E k )/E k = 0, k=l,...,N-l 

and considers the relation T{E)^\ = exp(P(E))^/i at the points Q k = (±, E k ) on the Riemann 
surface situated above E k 

<Pi{Qk)\ _ nP (E k ) ( <Pi(Qk 




xiiQu)) \xi(Q k )J- [ ] 

Solving it one obtains the values of P(E k ) on two sheets of T N . If the point Q(E k ) belongs to 
the same sheet of on which D(E k ) = e p( * Ek \ then one gets from ( |2.20| ) that pi(Q k ) = 0. 
Thus, the upper component of the Baker- Akhiezer function ^i(Q) has g + 1 zeros on above 
the points E k defined as roots of the polynomial B(E)/E. In similar way, g + 1 zeros of the 
lower component of ^i{Q) are related to the zeros of the polynomial C(E)/E defined by the 
C— component of the monodromy matrix. In general, two components of the Baker-Akhiezer 
function, ^"(Q; {t}) (a = 1,2), have different sets of g + 1 zeros on T N and their positions on 



Tjv depend both on the reggeon number, n, and the evolution times, {t}. Relation (|2.17|) implies 
that one of the roots should be at the points Qq. Indeed, it follows from fl2.14p and ( |2.7j ) that 



the polynomials B(E)/E and C(E)/E vanish at E — for L_ = and L + = 0, respectively, 
and therefore the components of the Baker-Akhiezer function vanish at the points Q$ = (±, 0) 
in agreement with (|2.13| ). 



The remarkable property of the roots E k = E k (zf\pf \ z^\p^) and the corresponding 
values of P(E k ) is that they form the set of separated coordinates for the system of N reggeons 
H|. Namely, they belong to the spectral curve 



e P(E k ) +e -P(E k ) = A( £ fe) = 2 + g2E 2 + + QnE N ; k = h ; N _ x (2 21) 

and have the following Poisson brackets 

{E k , E n } = {P(E h ), P(E n )} = , {E k \ P(E n )} = 5 kn . 

Inverting the relations fl2.21p one can express q k in terms of E k and then obtain the equations of 
motion for zeros E k of the Baker-Akhiezer function generated by hamiltonians q k in the following 
differential form M 

d ** = -£ / j { k = 2,3,...,N. (2.22) 

These equations can be integrated by the Abel map and their solution describe linear reggeon 
trajectories on the Jacobian variety of the Riemann surface [H]- As we will show in Sect. 3, 
using the Baker-Akhiezer function one can perform the inverse Abel transformation and construct 
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the explicit expressions for the reggeon trajectories on 2— dim plane of transverse coordinates in 
terms of Riemann theta-functions. 

For the Baker- Akhiezer function to be a well defined meromorphic function on the Riemann 
surface the number of its zeros, g + 1 = N — 1, should match the number of simple poles, which 
we choose to be at the points 71, 7 ff+ i in general position on r^. Moreover, considering 
the relations Q2.6|) and ( p. 13 ) at the points Q = (±,E) close to 7^ on Tn one finds that the 
components of ^%(Q; {£}) share the common set of poles and the positions of 71, 7^+1 on Tn 
do not depend neither on n, nor on times t. This set can be considered as part of the initial data 
for the evolution equations (|2.1|) . 

Finally, let us examine the relation ( |2.13| ) i n the neighborhood of punctures = (±, 00) on 
two sheets of Tn- Taking the limit E — > 00 one obtains from ( |2.13| ) that the solutions have the 
following asymptotics on, say, (— ) sheet, <p n (Q Q^o) = 0(E n ) and Xn(Q — > Q^o) = 0(E n ). 
Let us require that the solutions to ( |2.13|) should have a similar behaviour on the upper sheet, 
Vn{Q -> Qto) = 0(E~ n ) and Xn(Q -> Qto) = 0(E~ n ). This behaviour is consistent with (|2T3| ) 
provided that the r.h.s. of ( |2.13|) scales as E° at large E, or equivalently XniQto) ~ z nfn(Qto) = 
leading to the expression for the reggeon coordinates as a ratio of the components of the Baker- 
Akhiezer function at the point Q^. Thus, having constructed the Baker- Akhiezer function 
satisfying desired asymptotic behaviour one will be able to determine the reggeon coordinates. 



3. Construction of the Baker- Akhiezer function 

The existence of the Baker- Akhiezer function with the properties established in the previous 
section can be deduced from the following well-known fact in the theory of finite-gap soliton 
solutions |1(J. For a smooth hyperelliptic algebraic curve of genus g = N — 2 with 2 punctures 
at = (±, 00) and a given set of g+ 1 points (71,72, ...,7 9 +i) in general position there exists 
a unique function \I/"(Q; {t}) such that 

— isa meromorphic outside the punctures and has simple poles at the points 71, 7 9 +i; 

— \I/„ satisfies the normalization conditions ( |2 . 1 3| ) at the points Qq = (±, E = 0); 

— At the vicinity of two punctures on the upper and lower sheets the components ty® (a = 1, 2) 

have the following asymptotics 



*n(Q - M) = E* n e 



■±(n.{ T \) + 0[± 



(3.1) 



where 0^ is some independent function. 

The following comments are in order. In this definition, the Baker- Akhiezer function depends on 
the set of parameters {r} = r 1; ttv-2? which are different, in general, from the evolution times 
{t} entering (|2.1|). However, as a function of {r}, ^"(Q; {t}) satisfies the system of first-order 
linear matrix equations similar to (|2.6| ) and, as we will show in Sect. 4.1, the evolution times 
{t} corresponding to the hamiltonian flows are related to auxiliary parameters {r} by a linear 
transformation, t, +2 = J2k ^jk^k- 

Finally, substituting expression ( |3.1| ) into ( |2.13| ) and comparing the asymptotic behaviour of 
different sides of the relation (|2 . 1 3| ) for Q — > we obtain the following consistency conditions 

to = <f>t(n,{r}) = 02 (n + l,{r» 

n 0^(n,{r}) 0r(n + l,{r}) { ' } 
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and 



(") 



>t(n,{r}) 



1 0r(rc + l,{r}) 
z ( Z +1 <Pt(n+l,{T}) z^_ hn 4>T(n,{r}) 



(3.3) 



Together with ( |2. 16| ) these relations provide the solution to the hierarchy of the evolution equa- 
tions for reggeon coordinates and momenta, ( |2.1| ), in terms of the Baker- Akhiezer function. 



3.1. Basis of differentials 

To write down the expression for the Baker- Akhiezer function satisfying the above conditions 
one defines the canonical basis of cycles on as shown in Fig. 1 and constructs the basis of 
normalized differentials as follows JTTJ . 



The unique set of holomorphic differentials on Tn, or differentials of the 1st kind, dujk (k 
1 N — 2), is defined as 



N ^ 2 TT dEE^ _ N ^ TT dE E j 

duJk — > U k j — — — — > U k j 

U yW U JA 2 (E)-A 



(3.4) 



The coefficients Ukj are fixed by the normalization conditions 



dcJh = 27r<L 



jk i 



j,k = l,...,N-2 



and they can be calculated as an inverse to the following matrix 



A 2 {E)-A 



(3.5) 



(3.6) 



The unique set of meromorphic differentials on Tn of the 2nd kind, dQ U) (j > 1), with the j—th 
order pole at the punctures is defined as 



dVL 



U) 



dE 



A 2 {E) -4 



3Qn 



E^ 1 + jq^E^ 2 + ... + 0{E) 



The coefficients in front of the remaining powers of E in the numerator are fixed by the normal- 
ization conditions 



dtt U) = . 



dtt 



±(jEi~ 1 + 0(l/E 2 ))dE. 



(3.7) 



The unique dipole differentials on T^, or differentials of the 3rd kind, dQ^ and dQo, having 
simple poles at the points and Qq, respectively, and normalized by the conditions 



dttoo = <t> dVtQ = , d£l c 
are defined as 

d^co ■ 



E->oo dE 



dE 

)± E V 



1 dEA'(E) 



dE 



N 



A 2 (E)-A JA 2 (E)-A 



-q N E N ~ l + ... + 0(E) 



(3.9) 
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and 

HhdE N ^ 2 nT , dE E j 

dVt = —== +i]>( U o)j—r ( 3 - 10 ) 



To satisfy the normalization conditions ( |3.8|) the coefficients (U )j have to be equal to 

h N ~ 2 r rIF 

(Uo)k = --£ U jk I . (3.11) 



A 2 (£) -4 

with the matrix [7^ defined in Q3.4| ) and (|3.6| ). 

Let us also define the following (iV — 2)— dimensional vectors V, W {j) and 

V h = -il dn^ , = -i / cfft Ci) , A k {Q) = [ Q duo k , (3.12) 

Hk " J Pk ho 

where k = 1, JV — 2, 70 is an arbitrary reference point on Tn and integration in A(Q) goes 
along some path on between the points Q and 70. 



3.2. Baker- Akhiezer function 

The components of the Baker- Akhiezer function, ty°(Q; {r}) (a = 1,2), satisfying the conditions 
formulated in the beginning of Sect. 3 can be expressed in terms of the theta- function defined by 
the Riemann surface Tm (for definition see Appendix A) as follows |L4| 



U) 



WQ. W ) K{Q) ° {A{Q) + VU + E ^ W< ^ + Za)6{Zo) cxp n f Q dn +Tr-f 
n[Q ' 1 !) ~ K{Q a ) 9{A{Q) + Z a )9(Vn + Ef.f W™t, + Z ) GXP ^ k + T > L 

(3-13) 

Here, Q\ = Qq and Q2 = Qq are normalization points, the (N — 2)— dimensional constant vectors 
Z a and Z are given by 

Z a = Z - A(Q a ) , Z = A(Qi) + A(Q 2 ) - £ A( ln ) - K 

n=i 

with /C being the vector of Riemann constants. The function h a (Q) is defined as 

, (n) = 9(A(Q) + Z a )9(A(Q)-T a ) 
naW) 9{A{Q)-T + )9{A{Q)-T.)' 

together with the constants 

T + = 5>( 7s )+,4(7 3 ) + /C, T_ = Y^A(y s )+A( lg+1 ) + )C, T a = £ A{Q S ) + £ A{ ls ) +K . 

8=1 8=1 J = !> 2 S = l 

The following comments are in order. According to ( |3.12j ), the definition of the vectors A(Q a ) 
and A(j n ) implies a certain choice of the integration contours connecting the points Q a and 7„ 
with the reference point 70. The same contours should enter into the definition of the contour 

integrals of meromorphic differentials in the exponent of ( |3.13| ), J~ dfl u) = f® dfl U) — f® Q a dfl U) . 
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Only under this condition the Baker- Akhiezer function ( [3.13|) does not depend neither on the 



choice of the integration contours, nor on the reference point 70. Indeed, deforming the integration 
contour between, say, points Q and Q a as 7^ — > 7^ + J2j ^j a j + Y^j^jPji an d using the 
transformation properties ( |A.3| ) of the Riemann theta-function one gets that the variations of 
the exponent and the prefactor in (|3.13|) compensate each other, while the function h a (Q) stays 



invariant. 

Let us show that thus defined Baker- Akhiezer function ( |3.13| ) satisfies necessary normaliza- 



tion conditions. Considering the expression for h\{Q) and using the properties of zeros of the 
functions, (|A.4), we notice that two functions in the denominator of h\{Q) vanish at the 



points 71, 7 s _i, 7 9 and 71, 7 s _i, 7 9 +i, respectively, while the numerator vanishes at the points 
71, ...,7 s -i,Q2 phis additional g points coming from the first 9— function. In the expression for 
the Baker- Akhiezer function, ( |3.13|) , the latter cancels, however, against the same factor in the 



denominator of ( |3.13| ) and its zeros are replaced by zeros of the first function in the numerator 
of (|3.13| ). Therefore, the component ^\{Q] {r}) has simple poles at the g+ 1 points 71, ...,j g+ i 



and it vanishes at the point Q2 = Qo plus additional g points Q\, Q g . Being solutions of the 
equation 9(A(Q S ) + Vn + E^~i 2 W U) tj + Z x ) = they satisfy the relation (|OD , or 

Af-2 N-1 N-2 

£ A(Q.) + A(Q ) - £ A( ls ) = -Vn-J2 W U) r 3 . (3.14) 



One checks that for Q —>■ Q 1 = Qq different factors in (|3.13|) cancel against each other leaving 
us with ^niQo) = 1- As to asymptotic behaviour of ^(Q) for Q — > Q^, it is controlled by the 
exponent in (|3.13j ) and it is in agreement with fl3.1|) due to normalization conditions ( |3.7|) and (|3.8|) 
for the meromorphic differentials. The generalization of this analysis to the component ^(Q) 
is straightforward. We conclude that ( |3.13| ) gives a unique expression for the Baker-Akhiezer 
function corresponding to the N reggeon system. 



4. Explicit solutions 

We examine the asymptotics of the Baker-Akhiezer function (|3.13| ) in the vicinity of infinity 
and compare it with ( |3.1| ) to find the explicit expression for the function 0+. Substituting it into 
( |3.2|) and ( |3.3|) we obtain the solution to the hierarchy of the evolution equations for the reggeon 
coordinates and momenta in the following form 

= p V , °-\ T - + —^- exp liVon + iT Wtf V< (4.1) 

and 



Pn Zn ' n+l qN e(vn + EfJi 2 + z ) e(v(n + 1) + eS 2 + z + ) ■ [ ] 

Here, the constant p is given by a ratio of functions and its value can be made arbitrary using 
the symmetry ( |2.1S| ). The (N — 2)— dimensional constant vector (Z + )k is equal to 

Z + = A(Q+) + Z 1= A(Q+) - A(Q+) + Z = [ Q °° du + Z . 

JQ+ 
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The phase shift A and the oscillator frequencies are defined as 



° du k , V = -i ° dQ^ , Wtf } = -i ° dn u) , (4.3) 
Qo J Qo J Q 

where the normalized differentials are integrated along the path on Tjv connecting the points 
Qq and going through the reference point 70. Although 70 enters into the definition ( 3.12|) of 
the vector A(Q), neither the Baker- Akhiezer function ( |3.13|) , nor the reggeon coordinates (|4. 1|) 
and ( |4.2|) depend on its particular choice. Prefactor, q]( N , in (|4.2|) originates from the following 
expansion 

and it is consistent with the definition ( OJ) of the hamiltonian q^, that is Ilk=iP^ z n^n+i = 
<7jv. The expressions ( |4.1| ) and ( |4.2| ) involve iV — 1 free parameters - prefactor p and (N — 
2)— dimensional vector Z + (or equivalently Z ). 

Let us show that solutions ( |4.1| ) and ( [4. 2| ) satisfy the periodicity conditions ( |2.2| ). Using 
transformation properties of the 9— function, (|A.3| ), the same conditions can be expressed as 

V k = -i f dQ^ = 2tt— , V = -i dtl^ = 2vr— , 

with irtk and m being integers. To verify them we take into account the relation between eigen- 
values of the monodromy matrix, (|2.9|), and the dipole differential (|3.9|) 



dP(E) = —NdQoo , (4.4) 

and calculate from ( p.9|) and ( 2.11|) the multivalued function P(E) at the branching points Cj of 
the curve as e p ^ = ±1 and at the points Qq above E = as e p (°) = 1. 



4.1. Evolution times 



We recall that parameters r k entering (|4. Ij) and (|4.2|) do not coincide with the evolution times tk 
in (|2.1D. To find the relation between them we notice that g + 1 points Q\, Q g , Q g +i = Qq 
obey (|3.14|) being zeros of the component ^(Q) of the Baker- Akhiezer function. Moreover, 
the same points correspond to zeros of the B— component of the monodromy matrix and their 
t— evolution is described by ( |2.22| ). Comparing the relations ( [3.14 ) and Q2.22| ) and using the 
definition (ft.!2|) of the vector A{Q S ) together with ( p.4|) one obtains that up to a unessential 
additive constant 



N-2 



N-2 



E wjpTj = E u kj t 



i+2 j 



(4.5) 



where matrix U was defined in (O) and ( p.6|) . Let us also consider the r— dependent part of the 
exponent in (4.1) and take into account the relations ( |B.3|) and ( |4.5|) to get 



N~2 



zE w t 



T; 



3=1 



ih N ~ 2 

~ E wj" 

n 3,k=l 



k T j 



dE 



a k JA 2 (E) - 4 



7T 



V-2 

E U kj f 



dE 



j,k=i 



A 2 (E) 



tj+2 ■ 
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Comparing this identity with (|3.11 



) and ( |4.3|) one obtains 



AT-2 N-2 

E WPrj = E (Uo)jt j+ 2 ■ (4.6) 

3=1 3=1 



Substituting transition formulas ( flop and ( |4.6|) into ( |4.1|) and (|3.13| ) we restore the dependence 
of the reggeon coordinates and the Baker- Akhiezer function on the evolution times tk- Finally, 
combining together ( f4.1|) , ( [4.5| ), ( |4.6| ) and ( [2.19| ) we obtain the expressions for the holomorphic 
reggeon coordinates as 

z*\t) = pexp (2iht 2 + iV n + iU • t) ^ + U ' * + Z +7 P , (4-7) 

6(V n + U ■ t + Z + ) 

where the following notations were introduced 

N-2 N-2 
Uq ■ t = ^ (Uo)jtj+2 , U k -t= Ukjtj+2 ■ 

3=1 3=1 

The ratio of ^—functions in ( |4.7[ ) describes linear trajectories on the Jacobian JiT^) of the curve 
fl2.10|) . Relation (|4.7|) implies that the motion of reggeons becomes an image of these trajectories 
on the complex z— plane modulated by U(l) rotations with the basic frequencies (Uo)k- 

4.2. Curves with real branching points 

In the above consideration, the conserved charges qk as well as the branching points of the 
curve r^v were assumed to have arbitrary complex values. Let us consider the situation when all 
2N — 2 branching points (|2.11|) of the algebraic curve are distinct and real, 



ei < ... < e 2N -2- (4.8) 



The branching points were defined in ( 2.11|) as roots of the polynomial A (E) — 4 and in order 



for ek to be real the values of charges qk should also be real and, most importantly, they have to 
satisfy certain additional conditions. In particular, q 2 should be negative 

<? 2 < , (4.9) 

while the explicit form of the conditions on q 3 , q^ depends on the number of reggeons N 
[[]. p. For example, for iV = 3 states, they can be expressed as 

0<u 2 3 <^, (4.10) 

and for N = 4 states as 



_ 8 « 4 V48M4 + 1 - 2 
- 4« 4 < u\ < — \ -L- , (4.11) 

9 v48-u 4 + 1 — 1 

where w 3 = qzl[—q2) i ^ 2 and u 4 = q^jq\- Let us identify the moduli of the hyperelliptic curve 
Tn as (N — 2) -dimensional vector with the components Uk = qk/(—q2) h ^ 2 (k = 3,...,N). The 
relations ( |4.10| ) and ( |4.11|) define the regions on the moduli space, A4(T]\r), corresponding to 
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curves Tn with real branching points. At the boundary of these regions two of the branching 
points merge and the curve becomes singular. 

There is a special interest in considering soliton solutions for curves with real branching points. 
We recall, that the reggeon trajectories (4.1) appear as solutions of the evolution equations 



generated by the leading term So of the WKB expansion of the wave function of iV reggeon state 
( |1.6|) . A natural question arises whether the WKB expansion gives a meaningful approximation 
to solutions of the Schrodinger equation flLip . It turns out that for the quantum numbers qk 
satisfying ( |4.10|) and ( [4.11|) one can construct the exact solutions to ( |1 . 1|) within the Bethe 



Ansatz approach [T5| . A throughout analysis shows that the WKB expansion is in a 
good agreement with the exact solutions. Thus, similar to electron orbits in the Bohr model of 
hydrogen atom, the reggeon trajectory, corresponding to curves Tn with real branching points, 



describe the quasiclassical limit of the quantum states constructed in f|, [L5| . 

Let us consider the properties of solitons (|4.1|) under additional condition (|4.8|) . According 
to our definition of the canonical set of cycles shown in Fig. 1, the cycle aj lies on above the 
interval [e 2 j, e 2j+i]j called forbidden zones, on which A 2 (E) > 4. The Riemann surface defined 
by the curve r^v is constructed by gluing two copies of the complex E— plane along 2N — 2 finite 
intervals [e 2 j,e 2 j + i] and one infinite interval [e 2 N- 2 ,ei] going through infinity. On the intervals 
[e2j_i, e 2 j] (j = 1, N—l), called allowed zones, one has A 2 (E) — 4 < and the eigenvalue of the 
monodromy matrix, P(E), has pure imaginary values with ImP(E) being a quasimomentum. 
In particular, since A(0) = 2, the point E = belongs to one of the allowed zones, which we 
denote as [e 2K -i, &2k], 

ei < ... < e 2K -i < < e 2K < ••• < e 2N -2 • (4.12) 

The values of e 2 x-i and e 2 x can be defined as roots of ( |2.11| ) closest to the origin. The points Q$ 
belong to the different sheets of and one has to specify the integration path on Tjv between 
the points Q$ and Qq entering (|4.3|). We would like to stress that although its choice in (fO|) 
can be arbitrary, in order to apply the relations fl4.5| ) and ( f4.6|) one has to require that the path 
should not cross the canonical cycles aj and f3j. The later condition fixes it uniquely as shown 
in Fig. 1. The path is trapped between the cycles /3r-i and (3k and its orientation on T N is 
opposite to that of the (3— cycles. 

Under these definitions it is easy to see from (|3.6| ), and (ft.ll|) that parameters Ujk, (Uo)k, 



V, Vq and A are real. Moreover, according to ( [4.4D and (|3.12j ), the components of the vector V n 



are given by the difference of values of quasimomentum at the branching points e 2 Ar_2 and e 2n+ i, 
that is 



-%(N-n), l<n<K-l 
-^(N-n-1), K <n<N -2 



where integer K is defined by configuration of the branching points. In a similar way, the 
parameter Vq can be calculated as 

f e ^K r 2tT 27T , 

v = 2i dn^ + if dn 00 = — -v K = —(N-K). 

Jo J/3 K N N 

For the vector A n we get 

A n = 2i / du n + i f dlO n = An' + 2lliTKn 

Jo Jp K 

where tk u is an element of Riemann matrix, (|A.1| ), and An' = 2 Yl!j=i U n j J e 



V 4 - A2 (£)' 
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4.3. Solitons for N = 3 reggeon states 



Let us consider in more detail the properties of the soliton waves for N = 3 reggeon states with 
the quantum numbers q 2 = —h 2 and q 3 satisfying the conditions Q4.9| ) and ( 4.10| ). Choosing in 
( |4.10| ) the branch < q 3 /K 3 < ^7= one matches the branching points of ^ into ( |4.12| ) to get 
K = 1. The soliton solution (14. 71) takes the form 



Zn (t 2} £3) 



ln .. x 4i7r , TT \9(-fn + Ut 3 + Z + + iA) 
p exp ( 2 ih t 2 + — n + iUo t 3 ) e{ _^ n + Uh + z+) 



(4.13) 



where p and Z + are arbitrary complex parameters defined by initial conditions. The values of 
U , U and A can be calculated from (|3.6| ), ( p.ll|) and (|4.3| ) as elliptic integrals 



U 



7T 



re 3 dEE 
Je 2 v^V^i 



-2h 



re3 dE 



> e 2 

e 3 dEE 
e 2 V / A 3 ^4 



A = 2U 



dEE 



1 V^A 2 



(4.14) 



with A(E) = 2 - /i 2 £ 2 + g 3 £ 3 and A(e x ) = A(e 2 ) = A(e 3 ) = -2. The modular parameter r 
entering the definition of the Riemann 9— function takes pure imaginary values, Imr > 0, 



1 



e 4 dE E 



' re 3 dEE 
J e 2 V / A^4 



(4.15) 



For q% inside the region ( f4.10| ), the expression ( 4.13|) defines quasiperiodic reggeon trajectories 
on the (x, y)— plane of transverse coordinates, z = x + iy. An example of such trajectory corre- 
sponding to the quantum numbers q 2 = —1 and g 3 = 1/7 is shown in Fig. 2. Interesting property 




Figure 2: The trajectory of the reggeon on (x, y)— plane described by eq.( 4.13 ) for t 2 = and 
n = 1. The parameters have the following values: p — 1, h — 1 and Z + = 3 — i. Dot denotes the 
initial position of the reggeon at £3 = 0. 



of the trajectories (|4.13|) is that at any moment of time t 2 and £3 the reggeon coordinates satisfy 
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the condition 

Zl-Za Z 2 ~Z 3 Z3-Z1 (zt - z 2 ){z 2 - z 3 ){z 3 - zi) .h 3 

I 1 = ~1 77 77 7 = % — • (4.16) 

z 1 + z 2 z 2 + z 3 z 3 + z 1 [z 2 + z 3 ){z 2 + z 3 ){z 3 + Zi) q 3 

It follows from the expression (|1.4p, q 3 = Zi 2 z 23 z 3 ipip 2 p 3 , after one excludes the momenta using 
L 3 = —h and L + = L_ = 0. 

For q 3 = and q 3 = h 3 / \^27 the branching points of the curve T 3 merge leading to r = and 
r = ioo, respectively. The curve T 3 becomes singular and the properties of soliton solutions are 
changed. 

4.3.1. Singularity at q 3 = ^= 

In this case, the branching points are located at e\ = —y/3/h, e 2 = e 3 = 2y3/h, e 4 = 3y3/h 
and the a— cycle on T 3 is shrinking into a point. Applying relations (|4.14j ) and ([4.1 5|) one gets 



h 2 

U = -U = r , A = -21n2, r^ioo. (4.17) 

v3 

We notice that in the limit r — > ioo and Z + = O(t ) the ^—functions in (|4.13| ) are replaced by 1, 
since only one term with m = survives in the infinite sum ( |A.2| ). The corresponding solution, 
z [ n '{t 2 ,t 3 ) = pexp (2iht 2 + — ^7^3) > describes three reggeons located in the vertices of 
equilateral triangle and rotating around the origin as a whole. In order to avoid such rigid 
reggeon configurations, one has to adjust the constant Z + as 



Z + = nr + Z 00 + 0( 



T 



and use the asymptotic behaviour of the function ( |A.2|) 

6{u + tit; r) = 1 + e~ lu , asr-> ioo . 
Finally, the soliton solution for q 3 = h 3 /\/27 becomes 

4°) (t 2 ,t*)=P e 2W ^ n / + , $ = —n - - ^00 (4.18) 
n w> 6) h 4(1+ e-**) 3 ^3 

with p and being arbitrary constants. This relation defines periodic reggeon trajectories on 
the plane and example of such trajectory is shown in Fig. 3. 

4.3.2. Singularity at q 3 = 

In this case, e% = —e 2 = —2/h and two remaining branching points merge at infinity, e 3 = = 
0(1/ q 3 ). Integration in ft4.14|) and (|+|) yields 

U = —irh 2 , A = itti , r — ► . 

To analyze (|4.13 ) at small r one uses the duality property of the 9— function [11 



0(z;r) = (-irr^exp (~)o( -; --) , (4.19) 



A.1XT } \T T 
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Figure 3: The reggeon trajectories on (x,y)— plane described by eq. (|4.18|) for t 2 = 0. The 
parameters have the following values: p = 1, h — 1 and = 1 + 2. Dots denote initial positions 
of reggeons at t 3 — 0. 



which maps r — > limit into asymptotic behaviour at r — > —ioo. Taking the constant Z + in the 
form Z + = Zq + rZi with Zo and Z\ arbitrary one gets from (|4.19| ) and ( |4.13| ) that in the limit 
r — ► the solution ( |4.13|) becomes £3— independent z ( n (^2^3) = (— ) n pexp(2/izt 2 + |Z ). To get 
a nontrivial solution one has to adjust the constant Z + as 



7T 



- • tZ. + G(r 2 ) . 



The resulting expression for reggeon coordinates looks like 
4°W3) = -4°W3)=pe 2 ^, ^--w 



^"(-/•i./,) = : ; ri/2,0)tanh(yt 3 -^i). 



(4.20) 



Here, the positions of reggeons with the numbers n = 1 and n = 3 are frozen in time £3 and the 
n = 2nd reggeon moves between them, z^\t 2 , — 00) = ^f'*^, 0) and z^^, 00 ) = ^{h, 0). It 
is easy to see that the trajectory ( [4.20|) is consistent with (|4.16| ). 



5. Multi-soliton solutions 

As we have seen in the previous section, the soliton solutions corresponding to the degenerate 
singular spectral curve can be written for iV = 3 in terms of trigonometric functions. In the 
finite-gap theory they are called multi-soliton solutions 0. Let us show that the same property 
holds for any N. Although multi-soliton solutions can be obtained as limits of a general periodic 
solutions (|4~7| ) for degenerate curve T^, there is a simpler way to write them directly without 
going through a complicated degeneration procedure. 

The multi-soliton solutions correspond to a degenerate case, when all but two roots of the 
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equation ( [2.11| ) are double, that is 



N-2 

A 2 (E) -A = q 2 N {E- eJ)(E - e + )E 2 J] (E - e 3 ) 2 , (5.1) 

3=1 



with e_ 7^ e + . In this limit, the branching points of the curve merge and the definition ( |2.10| ) 
of Tjv takes the form 

N-2 

lt ns : y 2 = (E-e„)(E-e + ), y = y x q N [J (E - e 5 ) . (5.2) 

3=1 

Here, y is a rational curve and the Riemann surface corresponding to r^ ng has a genus g = 0, a 
sphere. It is constructed by gluing together two sheets along the cuts (oo, e_] and [e+, oo). One 
can introduce a global complex parameter k on and parameterize the solutions to (|5.2|) as 



£ = /3f/c+~) +a, y = /3( K -^\ (5.3) 

with /3 = ^(e + — e_) and a = |(e+ + e_). In this parameterization, the points k = oo and k = 
correspond to the punctures Q^, on the upper and lower sheets, respectively. The permutation 
of sheets of Tm corresponds to involution in k— plane, n — > 1/k, that does change the value of E, 

E(k) = E(1/k). (5.4) 

All meromorphic functions (and differentials) on Tjv introduced before can be now rewritten as 
functions of the parameter k only. In particular, consider the function w = e p ( E \ the Bloch 
multiplier, defined in (|2.<j|). According to ( |2.12| ), it has a pole of order N at and a zero of 
degree N at Q^. Therefore, as a function of k it has a unique form w = an N . At the vicinity of 
the punctures one substitutes n ~ (3/ E — > and k ~ E/j3 — > oo into w(k) and matches the 
asymptotics of w(E) with ( |2.12|) to get the values of the constants, a = 1 and (3 = q N , leading 
to 

w = e p W = K *. 

To find positions of the branching points of T^ ng we substitute this relation into ( |2.9| ), 

A(E) =k n + k~ n = 2 + q 2 E 2 (K) + ... + q N E N {K) , (5.5) 

and notice that, according to (|5.1|) , A 2 (E) = 4 at the branching points and A(E) = 2 for E = 0. 
Thus, k 2N = 1 at the branching points e± and ej, while k n = 1 at two points Qo on T^ 118 above 
= 0. These conditions define 2 A" points on the complex k— plane 

Kj = e&, (j = l,...,2iV-l) (5.6) 

and we identify the corresponding values E(kj) = 2q N x ^ N cos J + a as 2A" roots of A 2 — 4 in 
( |5.1|) . Since = K2N-j, it follows from (|5.4|) that, in agreement with (|5.1|) , 2N — 2 roots are 
degenerate, E(kj) = £?(k 2 jv-j) an d j 7^ , N. The remaining two roots, j = , N, correspond to 
the branching points e + = E(k ) and e_ = E(k n ). Among AT — 1 pairs Kj, n 2 N-j of the double 
roots one has to specify the pair that corresponds to two points Qq on r^ ng above E = 0. Let 



19 



3 — jo be the corresponding index. Since e± ^ E(Kj ) 
constrained to be positive even 



and k^ o 



0< Jo <N, j /2 = Z. 



1 the values of jo ar e 



(5.7) 



This index parameterizes different multi-soliton solutions. For given j one solves the equation 



E(k 
as 



30) 



to find the constant a 



-2q 



-l/N 



N 



cos(7r j /iV) and obtain the branching points of r^ ng 



COS 



7TJ 

N 



cos 



V N 



e± = 2q 



-l/N 



N 



±1 — cos 



V N 



with j = 1, N — 1 and j ^ jo- The branching points fix uniquely the curve r^ ng and allow 
to evaluate the corresponding values of the quantum numbers q^. Substituting the first relation 
( |5.3|) into ( |5.5|) and comparing the coefficients in front of different powers of k in the both sides 
of (|5.5[) one obtains the system of equations on q k and, in particular, 



N 2 



4sin 2 (^ 



N-l 
JV 



2iVcos 



/TTjo 

V N 

-2/3 



—3 and 



M 4 V2 



For N = 3 and iV = 4 the solution to ( |5.7| ) is j = 2 leading to q 2 q 3 
—4, g 3 = 0), respectively. 

Let us consider the Baker-Akhiezer function ^ on r^ ng . Its properties formulated in Sect. 3 
imply that being functions of complex k it can be written in the form 



^(K;{T}) = «-»e 



><?iv 



if R a (K;n,{r}) 
Ro(k) 



(5.1 



where R® and R a (a = 1,2) are polynomials of degree N — 1 in k. The zeros of Rq(k) become 
poles of and in what follows can be considered as parameters of the multi-soliton solutions, 
Ro(k) = YlfSi^n — 7 S ) . To verify the asymptotic behaviour ( |3.1|) one has to take into account 
the relation 



E{k) 



-l/N 



k H 2 cos 

K 



(5.9) 



and perform the limits k — > and k — > oo. The unknown so far polynomials R a are uniquely 

(5.10) 



defined by the normalization conditions ( 2.17 ), 



JO) 



and additional N — 2 linear relations 




1 



l,...,iV-l 



J 7^ Jo, 



(5.11) 



The meaning of the last condition is that parameters Kj and 1 / Kj define the same point on the 
curve. It belongs simultaneously to both sheets of r^ ng and the value of the functions should 
be the same. Let us choose the polynomial R2(k] n{r}) in the following form 



-R 2 (k; n, {t}) = A x (ac — 1/aCj, 



JO) 



N-l 

ik« 



JV-3 



/%) + (AC 



Jo , 



a i 

1=0 



K 



JV-3-Z 
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with A and a\ (I — 0, N — 3) being some functions of n and r. One verifies that ^ satisfies 
the normalization conditions (|5.10|) provided that A is given by 



A(n,{r}) 



0K^ j j 



30 ' K j) 



XK jo e 



const x exp 



m- 



N 



N-2 

2z 2J r s si 11 



\ N J 



Qn 



i/N 



The conditions ( |5.11|) lead to the linear system of equations for ai that can be expressed as 

nti 1 (i - w) 



N-3 

Yl °« sin 

s=0 



2 3+ N 



e2 3 



1 KjKj 0/ 



N-2. 



(5.12) 



Here, the parameters Kj are given by ( |5.6|) and the phases are defined as 



$ i (n,{r})=n 



2ttj 
N 



N-2 

4 ^ r s sin 

s=l 



7TJS 

AT 



q N + Zj 



(5.13) 



with Zj being some constants depending on zeros 71, 7^-1. Solution to Q5.12|) can be written 
explicitly as a ratio of determinants. Comparing fl3.1|) with (|5.8| ) in the limit k — > 00 we get 



02 ( n 5 { r D = ^4 x (1 + a ) = const x j o 1 + a (e 



- 1 ) 



(5.14) 



where ao is rational function of all phases except e l * J 'o . Repeating similar analysis for the poly- 
nomial R\ one arrives at the relation 



f(n, {r}) = A 1 x (1 + a ) = const x e 2- 



l + aofe-^V^e-**' 



- 1 ) 



(5.15) 



with a being a function complex conjugated to a . 

To finish description of the multi-soliton solutions one has to turn to the evolution times tk 
using the transition formulas (|4.5| ) and (|4.6|) . Let us consider the meromorphic differential of 
the second kind entering the definitions fl3.12| ) and ( f4.3p . It follows from ( |5.9| ) and ( p.7| ), that on 
the degenerate curve T s ^ g the differential dQ U) has a AT— th order pole at the points k = and 
k, = 00. Therefore it can be written as 



dftW = q N j/N d( K 



with the prefactor fixed by the asymptotics ( |3.7|) . Substituting this identity into ( |3.12j ) and 

one gets 



W, 



U) 



A -3/ N ■ 

4:q N sm 



irjk 



W^ 3) = 4g^' /iV sin 



V N J 



(5.16) 



where k — 1, N — 1 and k 7^ jo- To evaluate the r.h.s. of ( |4.5| ) and ( |4.6| ) one has to examine 
the E— expansion of the differentials duk and <if2 on r^ ng and identify the coefficients £/jy and 
(Uo)k- We observe that due to (|5.1| ) and Q3.4| ) the differential dujk has N — 2 additional poles 
on r^ ng located at the points E = ej. The normalization condition ( |3.5|) implies that it has a 
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vanishing residue at all points except of E — e^, where its residue is equal to i. This gives the 
following system of equations on Ukj 



N-2 

E UktE* = 2q N 

3=1 



1-1/N 



sm 



'nk 
~N 



N-l 



(5.17) 



valid for any E. The analysis of the differential cffi can be performed in a similar way with the 
only difference that, according to Q3.8Q , the residue of dfl vanishes at all points E = ej (j ^ j ) 
and its residue at E — ej = is equal to 1, 



N-2 , • x N-l 

2h + £ (U^E* = 2q^ N sin (-^) x J[{E - e,) . 



(5.1* 



i=l 
'/JO 



The relations ( |5.17| ) and ( |5.18| ) allow to calculate the coefficients C/jy and (t/o)fe in terms of the 
branching points and, in particular, 



N-2 



2 sin 



[7Tjo\ 1 1/iV 



fc,JV-2 



2 sin 



7T 



1/JV 



One checks, that for = 3, leading to jo = 2, these expressions are in agreement with (|4.17|) . 
Observing the relation Uj 0t k = (Uo)k, we substitute (|5.16|) into the transition formulas ( |4.5| ) and 
( f4.6|) to obtain the phases ( |5.13|) in the following form 



$k(n, {*}) =n 



2iik 



N-2 



J+2 



1< k < N - 1 



3=1 



with the coefficients Ukj calculated from (|5.17|) . Finally, substitution of ( ^.14[) and ( |5.15 ) into 
yields the multi-soliton solutions for the reggeon coordinates 



z { n 0) (t) =pe 



2iht2+i<b 



Ml 



) 



l + a (e- 4C15l ,...,e 



(5.19) 



with p being a constant. For iV = 3 this expression coincides with (|4.18 ). Solving ( 5.12 ) for 
arbitrary one obtains the multi-soliton solutions ( |5.19 ) in the form of rational functions of 
exponentials, e i. 



6. Conclusions 

In the present paper we studied the asymptotic solutions of the Schrodinger equation ( |1.1| ) for 
the color-singlet reggeon compound states in multi-color QCD. In the Regge limit, a nontrivial 
QCD dynamics affects only transversal reggeon degrees of freedom and the Af reggeon compound 
states look like a system of A^ interacting particles on the 2-dimensional plane of transverse 
reggeon coordinates. This system possesses a large enough family of conserved charges qk for the 
Schrodinger equation to be completely integrable. We identified the eigenvalue of the "lowest" 
charge q2 as a parameter of asymptotic expansion playing the role of the effective Planck constant. 

In the leading order of the asymptotic expansion, quantum fluctuations are frozen and quan- 
tum mechanical motion of reggeons is restricted to the classical trajectories driven by the "action" 
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function So- The conserved charges q k are replaced by classical functions on the phase space of 
N reggeons. They have a mutually vanishing Poisson brackets and generate the hierarchy of 
evolution equations for reggeon coordinates. For N = 2 reggeon states, solution to the evolution 
equations describes 2 reggeons rotating around the center-of-mass with the angular velocity given 
by the conformal weight h. For N > 3, integrals of motion g 3 , generate new modes of 

the classical reggeon motion which we identified as soliton waves propagating on the chain of 
TV particles with periodic boundary conditions. Applying the methods of the finite-gap theory 
we constructed the explicit form of the reggeon trajectories in terms of Riemann theta-functions 
and studied their properties. 

The orbits of the classical reggeon motion are parameterized by the eigenvalues of the con- 
served charges q k (k > 3). They define the moduli of the hyperelliptic curve and appear as 
parameters of the soliton solutions. In the leading order of the asymptotic expansion the eigen- 
values of qu could take arbitrary complex values. Quantization of q^, qN appears as a result of 
imposing the Bohr-Sommerfeld quantization conditions on the classical orbits of reggeon motion. 
It was shown || that for N = 3 reggeon compound states the results of the WKB expansion for 
eigenvalue of q^ are in a good agreement with the exact expressions T3] obtained by means of 



the algebraic Bethe ansatz |16| for integer positive h. Solutions to the Bohr-Sommerfeld condi- 
tions || give the quantized values of qu as functions of the conformal weight h and some additional 
set of integers {/} = In~2) parameterizing different families of curves on the moduli space, 

Ik = Ikih; {/}) {k = 3, N). These curves can be found as solutions of the Whitham equations 
H describing the adiabatic perturbation [O] of the reggeon soliton waves. 
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Appendix A 



Let us recall the definition of the Riemann theta-function and some of its properties |TTJ. We 
choose the basis of cycles otj, /3j on the Riemann surface with the canonical matrix of inter- 
sections as shown in Fig. 1. Then, we construct the basis of normalized holomorphic differentials 
and evaluate g x g Riemann matrix of their (3— periods 

1 f If 

doo k = r jk , — f du k = 5 jk . (A.l) 



27r/ ft . J 1 2tt 



The t— matrix is symmetric and it has a positive definite imaginary part. The basic g— dimen- 
sional vectors (ej) k = 2ii5j k and (bj) k = 2ixTj k generate the lattice C and allow to construct the 
torus J(Tn) = € g /£ called Jacobian of the Riemann surface Tn- 
One defines the (/—dimensional Riemann theta-function as 

9(u;t)= 22 exp (i7r(rm, m) + m)) (A. 2) 



Z 9 



where u = (ui, ...,u g ) is a complex vector on the Jacobian <J(T N ), m = (mi, ...,m 9 ) is a vector 
with integer components, (rm,m) = Y^%k=i T jk m j m k and (u,m) = J2 9 j=i u j m j- 
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The theta-function has the following periodicity properties. For any integer I 6 7L 9 and 
complex u G (D 9 one has 

6(u + 2ixl) = 6{u) , e(u + 2nrl) = exp (-in(rl,l) -i(u,l))6(u) . (A.3) 

The zeros of the theta-function satisfy the following condition. Let us fix a reference point 7 
on T N and define the vector A k (Q) = f^duk for any Q 6 r^v- Then, for an arbitrary complex 
Z = (Zi, Z g ) e C 9 the function — Z) either vanishes identically or it has exactly g 

zeros at the points Qi, Q g such that 

A(Q 1 ) + ... + A(Q g ) = Z-K (A.4) 

with K = K,(Tn,jo) being the vector of Riemann constants. 



Appendix B 

Let us establish useful relations between parameters and W defined in ( |4.3| ) and periods 
of normalized differentials on T^. Their derivation is based on the following identity (Stokes 
theorem) valid for any two unnormalized differentials du and du' on 

f2( a i b 'j - a 'j b j) = f du'(Q) f Q dw = - I duo{Q) f Q du' . (B.l) 

Here, aj = § a . du, bj = §g.du and similarly for du' are periods of the differentials, 70 is a point 
on Ttv in general position, C is an oriented closed contour on Fat which encircles all singularities 
of differentials. The contour C should not cross the cycles a s and (3 S (s = l,...,g). The same 
condition is imposed on the integration path between points 70 and Q in ( |B.1| ). 

Applying the identity ( |B.1|) to a pair of the normalized differentials du k and dQo defined in 
and ( p.iop one uses normalization condition on the l.h.s. calculates the r.h.s. by taking the 



residue at the poles Qq and gets the well-known relation JTT 



duk = — % f dQo . 

Qo ' hh 

Let us consider the following unnormalized dipole differential 

2ihdE 2ihdE 



dQ 



E v( E ) Ja 2 (e) - 4 



which appears as a leading term in the expansion (|3.10| ) of the normalized dipole differential dfV 
It has simple poles at the points Qo with the residue ±1, respectively, and its asymptotics at 
the vicinity of puncture on the upper sheet is 

d Q E ^ —dEE- N (l + 0(1/ E)) . (B.2) 

Applying the identity ( |B.1| ) to a pair of dQo and normalized meromorphic differential of the 
2nd kind, dQ u> , one takes into account normalization condition ( |3.8|) , calculates the r.h.s. by 
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taking residue at the simple poles Qq of the differential <if2 and j—th order poles of the 
meromorphic differential dQ U) and arrives at the following relation 



<f> d&of dQ 0) = -2ni f 

k = l a k JPk JQ 



Qo 1 
dVl U) + 2-2wi 



d 



3-1 



(j - 1)! Uz 



fi2 



2 =l/£;^0 



We notice that the derivative in the last term vanishes for j = 1, N — 2 due to 
gets the identity 



and one 



U) 



-o ,^/-s ih 



ih N ~ 2 r 



dE 



k= i J "k JA 2 (E) - 4 



l<J<N-2 



(B-3) 



with the vector W (J) defined in ( |3.12| ). 
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